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Abstract 



A review is given of some recently obtained results on analytic contractions of Lie al- 
gebras and Lie groups and their application to special function theory. The contractions 
considered are from 0(3) to E(2) and from 0(2,1) to E(2), or E(l,l). The analytic 
contractions provide relations between separable coordinate systems on various homoge- 
neous manifolds. They lead to asymptotic relations between basis functions and overlap 
functions for the representations of different groups. 
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1 Introduction 



Lie algebra contractions were introduced into physics by Inonii and Wigner [1] in 1953 as a 
mathematical expression of a philosophical idea, namely the correspondence principle. This 
principle tells us that whenever a new physical theory surplants an old one, there should exist 
a well defined limit in which the results of the old theory are recovered. More specifically Inonii 
and Wigner established a relation between the Lorentz group and the Galilei one in which the 
former goes over into the latter as the speed of light satisfies c — > oo. 

The theory of Lie algebra contractions (and deformations) has acquired a life of its own. 
It provides a framework in which large sets of Lie algebras can be embedded into families 
depending on parameters. All algebras in such a family have the same dimension, but they are 
not mutually isomorphic [2]. 

Two types of Lie algebra contractions exist in the literature. The first are standard Inonii- 
Wigner contractions [1, 3, 4]. They can be interpreted as singular changes of basis in a given 
Lie algebra L. Indeed, consider a basis {ei, . . . ,e n } of L and a transformation = f/j/ c (e)e/ c , 
where the matrix U realizing the transformation depends on some parameters e. For e ^ 
(i.e. some, or all of the components of e vanishing) the matrix U(e) is singular. In this limit 
the commutation relations of L change (continuously) into those of a different, nonisomorphic, 
Lie algebra 11 . 

More recently, "graded contractions" have been introduced [5, 6, 7]. They are more general 
than the Inonii- Wigner ones and can be obtained by introducing parameters modifying the 
structure constants of a Lie algebra Li in a manner respecting a certain grading and then 
taking limits when these parameters go to zero. 

It is well known that there exists an intimate relationship between the theory of special 
functions and Lie group theory, well presented in the books of Vilenkin [8], Talman [9] and 
Miller [10]. In fact all properties of large classes of special functions can be obtained from the 
representation theory of Lie groups, making use of the fact that the special functions occur as 
basis functions of irreducible representations, as matrix elements of transformation matrices, as 
Clebsch- Gordon coefficients, or in some other guise. Recently, the class of functions treatable by 
group theoretical and algebraic methods has been extended to the so called g-special functions 
that have been related to quantum groups [11, 12, 13]. 

One very fruitful application of Lie theory in this context is the algebraic approach to 
the separation of variables in partial differential equations [14, 15, 16, 17, 18, 19]. In this 
approach separable coordinate systems (for Laplace-Beltrami, Hamilton- Jacobi and other in- 
variant partial differential equations) are characterized by complete sets of commuting second 
order operators. These lie in the enveloping algebra of the Lie algebra of the isometry group, 
or in some cases of the conformal group, of the corresponding homogeneous space. 

A question that up to the last few years has received little attention in the literature, is 
that of connections between the separation of variables in different spaces, e.g. in homogeneous 
spaces of different Lie groups. In particular, it is of interest to study the behavior of separable 
coordinates, sets of commuting operators and the corresponding separated eigenfunctions under 
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deformations and contractions of the underlying Lie algebras. 

A recent series of papers [20, 21, 22, 23, 24, 25, 26, 27, 28, 29] has been devoted to this 
new aspect of the theory of Lie algebra and Lie group contractions: the relation between 
the separation of variables in spaces of nonzero constant curvature and in flat spaces. The 
curved spaces were realized as spheres S n ~ 0(n + l)/0(n), Lorentzian hyperboloids H n ~ 
0(n, l)/0(n), or 0(n, l)/0(n — 1, 1). The flat spaces where either Euclidean E n , or pseudo- 
Euclidean E(n — 1,1) ones. The curved and flat spaces were related by a contraction of their 
isometry groups and the corresponding isotropy groups of the origin. 

The essential point of these articles was the introduction of " analytic contractions" . The 
contraction parameter is R which is either the radius of the sphere S n , or the corresponding 
quantity x\ — x\ — ... — x\ = R 2 for the hyperboloid. The contractions are "analytic" because 
the parameter R figures not only in the structure constants of the original Lie algebra, but 
also in the coordinate systems, in the operators of the Lie algebra, in the invariant operators 
characterizing the coordinate systems, in the separated eigenfunctions and the eigenvalues. 

Once the parametrization displaying the parameter R is established, it is possible to follow 
the contraction procedure R — > oo explicitly for all quantities. For the Lie algebra realized 
by vector fields, the Laplace -Beltrami operators, the second order operators in the enveloping 
algebras, characterizing separable systems of coordinates, the separated ordinary differential 
equations, the eigenfunctions and the coefficients of the interbases expansions. 

For two-dimensional spaces all types of coordinates where considered, for example, con- 
tractions of 0(3) to £'(2) relate elliptic coordinates on S2 to elliptic and parabolic coordinates 
on E2. They also relate spherical coordinates on S2 to polar and Cartesian coordinates on E2 
([20, 22, 24]). Similarly, all 9 coordinate systems on the H 2 hyperboloid can be contracted to 
at least one of the four systems on E 2) or one of the 10 separable systems on [21, 23]. 

Contractions from S n to E n where considered for subgroup type coordinates in Ref. [25, 26, 
28], for subgroup type coordinates and certain types of elliptic and parabolic ones. 

The main application of analytic contractions in this context is to derive special function 
identities, specially asymptotic formulas. Among other possible applications we mention the 
theory of finite dimensional integrable and superintegrable systems [30, 31] . 

In this paper we restrict ourselves to two-dimensional spaces of constant curvature. The 
organization of the article corresponds to the table of contents. 

2 Separation of variables in two dimensional spaces of 
constant curvature 

2.1 Operator approach to the separation of variables 

Let us first consider a quite general two dimensional Riemannian, pseudo-Riemannian or com- 
plex Riemannian spaces with a metric 





(2.1) 
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In this space we introduce a classical free Hamiltonian 

H = gik(u)PiP k , (2.2) 

where pi = —dH/dtii are the momenta classically conjugate to the coordinates u l . We also 
introduce the corresponding Laplace-Beltrami operator 

A ^ / — ^ ^ (*2 3) 

We will be interested in two related question: 

1. What are the quadratic polynomials on phase space 

Q = a tk (u)p tPk (2.4) 
that Poisson commute with the Hamiltonian 

In other words, when do quadratic (in the momenta) integrals of motion exist? Respectively 
what are the second order Hermitian operators 

Q = {a l \u)d Ui d Uk } (2.6) 

(where the bracket denotes symmetrization) that Lie commute with the Laplace-Beltrami op- 
erator 

[Q, H] = QH - HQ = 0? (2.7) 

2. Do the Hamilton- Jacob i and Laplace-Beltrami equations in the considered space allow 
the separation of variables, and if so, how do we classify and construct separable coordinates? 
By separation of variables for the Hamilton- Jacobi equation we mean additive separation 

dS dS 

^t-t- = A (2.8) 

S = S 1 (u l ,\,fi) + S 2 (u 2 ,\,fi). (2.9) 

For the Laplace-Beltrami operator we have in mind multiplicative separation 

A^ = A^, ^ = ^i^A^^^A,^). (2.10) 

In both case A and [i are the separation constants. 

In this review article we shall mainly be interested in Laplace-Beltrami operators in differ- 
ent spaces. However, some aspects of separation are simpler to discuss for the Hamilton- Jacobi 
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equation. In two-dimensional Riemannian space the two equations separate in the same coor- 
dinate systems. 

The existence of integrals of motion that are either linear or quadratic in the momenta was 
analyzed by Darboux [32] and Koenigs [33] in a note published in Volume 4 of Darboux's lec- 
tures. In particular it was shown that a metric (2.1) can allow 0, 1, 2, 3 or 5 linearly independent 
quadratic integrals. The case of 5 quadratic integrals occurs if and only if the metric corre- 
sponds to a space of constant curvature. In this case the second order integrals are "reducible". 
That means that the metric allows precisely three first order integrals 

Li = a i (u)p 1 + bi(u)p 2 , u=(u 1 ,u 2 ), i = 1,2,3 (2.11) 

and all second order integrals are expressed as second order polynomials (with constant coeffi- 
cients) in terms of the first order ones: 

3 

Q — A ik LiL k A ik = constant (2.12) 

i,k=l 

If the polynomial Q is the square of a first order operator L, then it will provide a subgroup 
type coordinate. This is best seen by considering the corresponding first order operator 

X = C(ui,u 2 )d Ul +r)(u 1 ,u 2 )d Ul (2.13) 

that generates a one dimensional subgroup of the isometry group G. /,From (ui,u 2 ) we trans- 
form to the new coordinates (vi,v 2 ) "straightening out" the vector field (2.13) to the form 

X = d Vl (2.14) 

Then v\ will be an ignorable variable. The complementary variable v 2 = <f)(ui,u 2 ) can be 
replaced by an arbitrary function of v 2 , the ignorable variable v\ can be replaced by f(vi)+g(v 2 ) 
where both / and g are arbitrary. The separable coordinates are V\ and v 2 (with the above 
mentioned arbitrariness). 

Now let us assume that an irreducible quadratic integral Q as in (2.4) is known for a 
considered metric (2.1) (that is, Q is not square of a linear integral). We can then impose that 
two equations be satisfied simultaneously. In the classical case they are 

OS dS x dS dS 

9ik^-T- = A, a ik — — = n. (2.15) 

Similarly, we can consider the quantum mechanics of a free particle in such a space and write 
two simultaneous equation 
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Separable coordinates for the two systems (2.15) and (2.16) are obtained by simultaneously 
transforming H and Q to a standard form, in which the matrices and of (2.15) (and 
(2.16)) are diagonal. This can be done by solving the characteristic equation 

kfc ~P9ik\ = 0- (2.17) 

If two distinct roots p 1 and p 2 exist, they will provide separable coordinates, at least over the 
field of complex number. If we are considering real spaces, then it may happen that p 1 and p 2 
are real only in part of the space and do not parametrize the entire space. We will see below 
that this indeed happens for instance in the pseudo-euclidean plane E\j_. 

The roots p\ and p 2 can be replaced by any functions u = u(pi), v = v(p 2 ). This freedom 
can be used to transform H and Q simultaneously to the form 

H = „M 1 W, \ (Pv + Pv) = A > Q= ^ 7 \ W u )pI - a ( v )pD = ( 2 - 18 ) 
a(y) + p{u) a[v) + p[u) 

The Hamiltonian H in (2.18) is in its Liouville form [34]. The separated equations are 

aH + Q = a\ + p, (3H - Q = (5\ - p (2.19) 

for the Hamilton- Jacobi equation and similarly 

(aH + Q)^ = (a\ + n)V, ((3H - Q)V = ((3\ - p)V (2.20) 

for the Laplace-Beltrami equation. 

Let us now restrict ourselves to two-dimensional spaces M of constant curvature, that is 
to the Euclidean plane E 2 , pseudo-Euclidean plane £1,1, sphere ^2 and two, or one- sheeted 
hyperboloid H 2 . Each of these has a three-dimensional isometry group G. The Lie algebra L 
of G has in each case a standard basis which we denote {X±, X 2 , X 3 }. 

The Laplace-Beltrami operator H = Alb (2.3) is in each case proportional to the Casimir 
operator C of the Lie algebra L. The operator Q commuting with H will have the form 

3 

Q = X! A ik XiX k A ik = A ki (2-21) 

i,k=l 

where A is constant matrix. Let g G G be an element of the Lie algebra of the isometry group 
of the considered space. Let us rewrite eq. (2.21) in matrix form 

Q = u T Au, u T = {X u X 2 , X 3 } (2.22) 

The transformation g acting on the space M induces a transformation u' = gu on the Lie 
algebra L. The Casimir operator Alb stays invariant, but Q transforms to 

Q = u T g T Agu. (2.23) 
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Thus, for spaces of constant curvature a classification of operators Q commuting with H reduces 
to a classification of symmetric matrices A = A T into equivalence classes under the congruence 
transformation 

A' = g T Ag, g e G. (2.24) 

This problem, as we shall see below, can be reduced to that of classifying elements of Jordan 
algebras into equivalence classes. 

Furthermore, the operator C of L can also be written in the form 



/ ci \ 

C = u T Cu, C= c 2 

V C 3 / 



The matrix A of equation (2.22) is 

a b T \ 

b s r ' \ c 3 c 2 



and C of (2.25) is 



(2.25) 



Two matrices A and A will give equivalent coordinate systems, if they satisfy 

A = \g T Ag + fiC, A ^ 0, (2.26) 
where A and ji are real constants. 

2.2 Separable Coordinate Systems in the Euclidean Plane. 

The Lie algebra of the isometry group E(2) is given by 

L = x 2 d Xl - x x d X2 , Pi = d Xl , P 2 = d X2 . (2.27) 
The operator Q of eq.(2.21) will in this case be 

Q = aL 2 + 6i(LPi + PiL) + b 2 (LP 2 + P 2 L) + Cl P? + c 2 P| + 2c 3 PiP 2 . (2.28) 
An P(2) transformation matrix will be written as 

0=(J r)> e T = (6,6), RelR 2 , R T R = I (2.29) 



A=(t bT Q ), S= ( J J V b T = ( bl ,b 2 ), (2.30) 



C= | 1 | , (2.31) 



since the Casimir operator is 

C = A = Pi + P 2 . (2.32) 

The transformation (2.26) with A = 1 has two invariants in the space of symmetric matrices A, 
namely 

h = a, h = {[a(ci - c 2 ) - (bl - bl)} 2 + 4(ac 3 - &i6 2 ) 2 }* (2.33) 
Correspondingly, the operator Q can be transformed into one of four canonical forms 

1) h = 0, I 2 = Q c = Pi, (2.34) 

2) h ^ 0, J 2 = Qfl = Z 2 ; (2.35) 

3) h = 0, 7 2 ^ g P = LP 2 + P 2 L; (2.36) 

4) A ^ 0, J 2 ^ g E = L 2 - D 2 P 2 , D 2 = ^ (2.37) 



The first two forms correspond to subgroup type coordinates. Thus, Q c of (2.34) corresponds 
to Cartesian coordinates, in which Pi = d x (and also P 2 = d y ) is already straightened out. Both 
x and y are ignorable. The second, Qr corresponds to polar coordinates 

x = rcos0 , y = rsm(p (2.38) 

in which L = is straightened out so that is an ignorable variable. 

The coordinates corresponding Q P of (2.36) are the parabolic coordinates 

x = \{u 2 -v 2 ), y = uv. (2.39) 

Equivalently, if we take Qp = LP\ + P\L the prescription (2.17) leads to 

x = uv, y = h u 2 -v 2 ). (2.40) 

Finally, Q E of eq. (2.37) corresponds to elliptic coordinates in the plane. They can be written 

as 

x = D cosh £ cos rj, y = D sinh £ sin rj. (2.41) 



2.3 Separable Coordinate Systems in the Pseudo-Euclidean Plane 

The Lie algebra of the isometry group E(l, 1) can be represented by 

K=(td x + xd t ), P = d t , P 1 = d x . (2.42) 
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Table 1: Orthogonal coordinate systems on two-dimensional Euclidean plane E 2 



Coordiante systems 


Integrals of motion 


oolution ol 










TT 1 "I 1 l_ J_ * 

Helmholtz equation 


1. bartesian x, y 








Exponential functions 


— oo < x, y < oo 


n — 
Use — 


p2 


p2 




TT T 4 . . 1 

11. rolar 








Product of 


x = r cos ip, y = r sin (p 


Qr = 


L 2 




Bessel function 


<r < oo, < if <2n 








and exponential 


III. Parabolic 








Product of two 


2 2 
ii — u 

x = 2 , y = uv 


Qp = 


LP 2 - 




parabolic cylinder 


< u < oo, —oo < V < oo 








functions 


IV. Elliptic 








Product of periodic 


x = D cosh £ cos 1], 


Qe = 


L 2 — 


P 2 P| 


and nonperiodic Mathieu 


y = D sinh £ sin 








functions 


< £ < oo, < ?7 < 2?r 











The second order operator (2.21) is 

Q = aK 2 + b (KP + P K) + b 1 (KP 1 + P 1 K) + c P 2 + Cl P 2 + 2c 2 P Pi. (2.43) 
Equivalently, the matrix A of (2.24) is 

A=(l b l), C=(H), e-frM (2.44) 

We will classify the operators Q into conjugate classes and the action of the group £'(1,1), 
including the reflactions 

n : (x,t) -> (x, -t), IT: (x,t) -> (-x,t). (2.45) 

An element of E(\, 1), acting on the Lie algebra (K, Pq, Pi) can be represented as 



T 



A 



with 



e = (£o,£i), A G IR , A JA = J (2.46) 
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Table 2: Orthogonal coordinate systems on the two-dimensional pseudo- Euclidean plane E^i 



Coordinate system 


Integrals of motion 


Solution of 
Helmholtz equation 


I. Cartesian 

t, X 

* 


Qc = PoA 


Product of exponentials 


II. Pseudo-polar 

t — r cosh Tii x = r sinh t 2 

r > 0, — oo < r 2 < oo 


Qs = K 2 


Product of Bessel 
function and exponential 


III. Parabolic of type I 
t = \{u 2 + v 2 ), x = uv 
V > 0, — oo < u < oo 


Q I P = {P 1 ,K} 


Product of parabolic 
cylinder function 
for t 2 - x 2 > 


IV. Parabolic of type II 
t = uv, x = \{u 2 + v 2 ), 
v > 0, — oo < it < oo 


Q$ = {P ,K} 


Product of parabolic 
cylinder function 
for x 2 — t 2 > 


V. Parabolic of type III 

t = Uv-C?-(v + 0, 
x = ^(v-C) 2 + (v + C), 

— OO < 7], ( < OO 


Qp 1 = { p o, K} + {Pi, K} 
+ (Po ~ Pi) 2 


Products of two linear 
combinations of Airy 
functions for x + t > 


VI. Hyperbolic of type I 
t = \ (cosh 2^ + sinh 2±C) 
x = | (cosh 2^ - sinh 2±C) 
— oo < rj, ( < oo 


Qh = K 2 — l 2 P P! 


Product of Mathieu equation 
solutions with argument 
displaced by in/ 2 


VII. Hyperbolic of type II 
t = (smh(ri - C) + e"+ c ) 
x = (sinh(?7 - C) - e^) 

— OO < Tj, ( < OO 


Q» = K 2 + {P 1 + P 2 ) 2 


Product of two solutions 
of Bessel's equation, 
one with real and one 
with imaginary arguments 


VIII. Hyperbolic Of type III 
t = (cosh(7] - C) + e T,+c ) 
x = (cosh(r] — () — e v+ ^ 
— oo < 1], ( < oo 


Qh 1 = K 2 - (Pi + P 2 ) 2 


Product of two solutions 
of Bessel's equation 
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Table 2: (Continue) 



Coordinate system 


Integrals of motion 


Solution of 
Helmholtz equation 


IX. Elliptic of type I 

t = D sinh i] cosh ( , x = D cosh r\ sinh ( 

— OO < Tj, ( < OO 


Q I E = K* + D 2 P? 


Product of two solutions 
of the nonperiodic Mathieu 
equation 


X. Elliptic of type II 

(i) t — d cosh rj cosh (, x = d sinh 77 sinh £ 

— 00 < 77 < 00, C > 

(ii) t — d cos i] cos C, x = <i sin 77 sin ( 
<r] <2tt,0 < ( <n 


Q» = K 2 — d 2 P 2 


Product of two solutions 

(i) of the nonperiodic 
Mathieu equation 

(ii) of the periodic 
Mathieu equation 



The matrix A of (2.44) is subject to the transformation (2.26) and in this case we have 

A'-n T An- ( ° a^ T + P T A \ (!? , R) 

A ~ 9 Ag ~ \ ia + A T [3 A T CA + A T ^ T + £,P T A + a^ T ) ' [ZA * } 

One of the invariant of this transformation is the constant a which can be closen to be a = 1, 
or is already a = 0. 

Let us first consider a 7^ 0. Choosing £ = — j3 T A and putting a = 1, we obtain 

A'=(lc,y C' = JA-V(C-/3/3 T )A. (2.49) 

Notice that C and C are symmetric matrices, but we have 

X = J(C - (3(3 T ), JX T = XJ, (2.50) 

that is, X is an element of the Jordan algebra jo{\, 1). Since A is an element of the Lie group 
0(1, 1), we are faced with a well known problem: the classification of elements of a Jordan 
algebra with respect to conjugation under the corresponding Lie group. The results are known 
for all classical Lie and Jordan algebras [35], and for jo(l, 1) they are quite simple. The matrix 
X can be transformed into one of the following 

with p G IR, q e 1R. 
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For a = 0, (&o,&i) 7^ (0,0) transformation (2.26) leads to eq. (2.48) in which we set a = 0. 
We choose the matrix A contained in 0(1, 1) to transform A T [3 to standard form and then 
choose £ to simplify the matrix C. For |6 | > |&i|, |&o| < |&i| an d |&o| — |&i| we can transform 
A into 







0\ 



2 -6? 








A, = 



V ^ - & o 

/ 1 1 \ 

1 7 _ 7 

V 1 -7 7 / 












o J 




bl) 










) 



±1 



±1 



7 = 0,1 



(2.52) 



(2.53) 



(2.54) 



respectively. 

Finally, for a = b = bi = 0, C ^ we can use A to transform C into one of its standard 
forms JXi, i = 1, 2, 3 with as in (2.51). 

Thus, we have obtained a classification of matrices A that determine the operator Q. Let 
us now list the corresponding operators. We first notice that if a = b\ = b<i = 0. The 
corresponding operator Q is in the enveloping algebra of a maximal Abelian subalgebra of 
e(l,l), namely (P ,Pi). Similarly, for a = 1 and X = X 1 in eq. (2.51) with p = q we find 
that Q = K 2 is in the enveloping algebra of a different maximal abelian subalgebra of e(l, 1), 
namely o(l, 1) (generated by K). These two cases correspond to subgroup type coordinates. 
The other ones to nonsubgroup type. 

The list of operators must be further simplifield by linear combinations with the Casimir 
operator 



C 



p2 



Pi 



(2.55) 



Finally, we obtain a representative list of 11 second order operators in the enveloping algebra 
of the Lie algebra e(l, 1). 

QxM) = a(P 2 + P?) + 2bP P 1 ; 

(a,6) = (l,0),(l,l),or(0,l) 
Q 2 = K 2 - 
Q 3 = KP. + P.K; 
Q 4 = KP + P K; 
Q 5 = K(P + P 1 ) + (P + P 1 )K; 
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g 6 = K(P + Pi) + (P + Pi)K + (P -P) 2 ; (2.56) 

q 7 = jr 2 -/ 2 p p i; i>o 

Qs = K 2 - D 2 P 2 ; D > 
Q 9 = K 2 -d 2 P 2 ; d>0 

Qio = K 2 + (P + P) 2 ; 

g n = k 2 -(p + Pi) 2 . 

To obtain separable coordinates we proceed as in Section 2.1. 

1. The operator Qi(a,b) for any a and 6, corresponds to Cartesian coordinates (t,x), since 
the operators that are really diagonalized are Pq and P\ (they correspond to a maximal Abelian 
subalgebra {P ,Pi} G e(l, 1). 

2. The operator Q2 = K 2 also corresponds to subgroup type coordinates, namely pseu- 
dopolar coordinates 

t = rcosho;, x = rsinhQ; (2-57) 

< r < 00, 00 < a < 00. 

These coordinates only cover part of the pseudo-Euclidean plane, since we have t 2 —x 2 = r 2 > 0. 
By interchanging t and x in (2.57) we can parametrize the part with t 2 — x 2 = — r 2 

The operators Q3, Qn can lead to separable coordinates via the algorithm of eq. (2.17). 

Two problems can and do occur. The first is that the roots of eq. (2.17) may coincide: p\ = p 2 - 
Then we do not obtain separable coordiantes. This happens in precisely one case, namely that 
of the operator Q5. 

To other problem that may occur is that the eigenvalues p(t, x) may be complex in at least 
part of the (x, t) plane. This part of the plane will then not be covered by the corresponding 
coordinates (px,p 2 )- 

The results of this analysis are presented in Table 2 and essentially agree with those of 
Kalnins [36]. 

2.4 The systems of coordinates on S2 

The Lie algebra of isometry group 0(3) is given by 

d 

Li = -e ikj u k —, [Li,L k ] = eikjLj, 1, k,j — 1, 2, 3. (2.58) 

where Ui are the Cartesian coordinates in the ambient Euclidean space E 3 . On the sphere S2 
we have u\ + u\ + u\ = R 2 . The Casimir operator is 

C = R 2 A LB = Ll + L 2 2 + L 2 3 (2.59) 

and the Laplace-Beltrami equation (2.10) for S% has the form 

A LB ^ = "^tJ^ P) = S, fc (a)$i fc (/3), (2.60) 
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where i = 0, 1, 2, .... The second order operator Q of eq. (2.21) is given by 

Q = AikLiLk, = Aki. 

The transformation matrix for 0(3) can be represented as 

(cos a cos (5 cos 7 — sin a sin 7 — cos a cos j3 sin 7 — sin ct cos 7 cos ct sin /3 
sin a cos /3 cos 7 + cos a sin 7 — sin a cos /? sin 7 + cos a cos 7 sin a sin /? 
— sin (3 cos 7 sin [3 sin 7 cos (3 

where (a, /3, 7) are the Euler angles. 

The matrix A ik can be diagonalized to give 



(2.61) 



(2.62) 



Q(ai, a 2 , a 3 ) = Q = a x L\ + a 2 L\ + a 3 Lg. 



(2.63) 



For ai = a 2 = 03 we have Q ~ 0. If two eigenvalues of are equal, e.g. a x = 02 7^ 03, 
or ai 7^ a 2 = 03, or ai = 03 7^ a 2 we can transform the operator Q into the operators: 
Q(0,0,1) = L 2 , (5(1,0,0) = L\ or (5(0,1,0) = L| respectively. The corresponding separable 
coordinates on S 2 are the three types of spherical ones 



Ml 

u 2 
us 



R sin 9 cos ip 
R sin 9 sin ip 
RcosO 



Rcos9' 

R sin 9' cos <// 

i? sin 9' sin 99' 



i?sin 9" sin <//', 

iZcosfl", 

i? sin 9" cos <//' 



(2.64) 



where p G [0, 27r), 9 G [0, 7r]. They correspond to the group reduction 0(3) D 0(2) and X = L\ 
is invariant under 0(2) and under reflections in all coordinate planes. 

The 0(3) unitary irreducible representation matrix elements of (2.62) results in the well- 
known transformation formula for spherical functions Yi m (9, ip) [8, 37] 



nm'(*V)= E D l mm ,{a,M)YU9, V ), 

m=—l 



(2.65) 



where Di 



nil m2 (a,(3, 7) - are the Wigner D-functions 



e~ ima d 



m,m 
m—m 



im 7 



(2.66) 



7 7TT\ 77 7T77 ;d( cos ^) SUl^f m 

(m-m')! \ (£-m)! (£ + m')! 2 2 



x F 



m — 1, —m — 
m — w! + 1 



; — tan 2 |/3 



(2.67) 



and the spherical angles (9, (p) and (9', ip') are related by 

cos#' = cos # cos /3 + sin sin /3 cos(</? — a) 

cot 6* sin (3 



cot(</?' + 7) = cot(y? — a) cos/5 



sin(y? — a) 



(2.68) 
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In particular, Yi m (6, (p) corresponding to the solution of Laplace-Beltrami equation in the sys- 
tems of coordinates (2.64) are related by the formulas 

1W*V) = E D l mm/ (0,~)YU0,<P), (2.69) 

YiM^^ = E Dl mm"(l^ Q ) Y l,m(0,(p), (2.70) 
m=—l 

= e ^"(o,f,f)n^V). (2.7i) 

m'=—l 

When all three eigenvalues aj are different, then the separable coordinates in eq. (2.60) are 
elliptic ones [38, 39, 40]. These can be written in algebraic form, as 

„,2 P 2 (Pi -Ql)(P2 -Ol) , p2 (Pl-a 2 )(p2-fl2) ,7 D 2 (Pi ~ <*3) (P2 ~ Og) (t) ^ 

u 1 = R - r, u 2 = R - -, u 3 = R - {2.. 12.) 

{a 2 - a 1 ){a 3 - aij [a 3 - a 2 ){a 1 - a 2 ) (ai - a 3 )[a 2 - a 3 ) 

with ai < pi < a 2 < p 2 < a 3 . 
In trigonometric form we put 

Pi = ai + (a 2 - ai) cos 2 0, p 2 = a 3 - (a 3 - a 2 ) cos 2 9, (2.73) 

and obtain 



ui = RV 1 — k' 2 cos 2 9 cos 0, w 2 = -R sin 6* sin 0, tt 3 = i?y 1 — k 2 cos 2 cos 9. (2.74) 

< < 2tt, < 9 < vr, 

where 

A; 2 = = sin 2 /, k' 2 = = cos 2 /. A; 2 + A;' 2 = 1. (2.75) 

a 3 - ai a 3 - ai 

The Jacobi elliptic version of elliptic coordinates is obtained by putting 

Pi = ai + (a 2 - ai)sn 2 (a, k), p 2 = a 2 + (a 3 - a 2 )cn 2 (/3, fc'), (2.76) 

We obtain 

■Ui = i?sn(o;, k)da(P, k'), u 2 — Rcn(a, k)cn((3, k'), u 3 — Rdn(a, k)sa(/3, k'), (2.77) 

-K <a<K, -IK' <(3< 2K', 

where sn(a, k), cn(a, k') and dn(/3, k) are the Jacobi elliptic functions with modulus k and k', 
and K and X' are the complete elliptic integrals [41]. 

The interfocal distance for the ellipses on the upper hemisphere is equal to 2fR. 
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Table 3: Orthogonal coordinate systems on two-dimensional sphere S 2 



Coordinate systems 


Integrals 
of motion 


Solution of 
Helmholtz equation 


Limiting 

systems on E 2 


I. Spherical 
ui = Rsin9 cosy? 
U2 = R sin 6 sin ip 
u 3 = RcosO 

e [0,2%), e e [o,tt] 


f~\ T 2 

= -^3 


Product of associated 
Legendre polynomials 
and exponential 


Polar 
Cartesian 


II. Elliptic 

Ui = Rsn(a, k)dn(/3, k') 
u 2 = Rcn(a, k)cn((3, k') 
u 3 = Rdn(a, k)sn((3, k') 
a E [-K, K], 
(3 e [-2K', 2K'\ 


Q E = k' 2 L 2 - k 2 L\ 


Product of two 
Lame polynomials 


Elliptic 
Polar 
Cartesian 
Parabolic* 



* After rotation. 

2.5 Systems of coordinates on H2 

The isometry group for the hyperboloid H 2 : u\ — u\ — u\ = R 2 , where Ui(i = 0,1,2) are 
the Cartesian coordinates in the ambient space £2,1 is 0(2,1). We choose a standard basis 
Ki, K 2 , L 3 for the Lie algebra o(2,l): 

Ki = - (u d U2 + u 2 d U0 ), K 2 = - (u d Ul + Ui<9 U0 ), L 3 = - (u^ - u 2 d ui ) 

with commutation relations 

[K 1 ,K 2 ] = -L 3 , [L Z ,K 1 \ = K 2 , [K 2 ,L 3 } = K 1 . 

The Casimir operator is 



C = R 2 A LB =K\ + K\-L\. 
and the Laplace-Beltrami equation (2.10) given by 

1(1 + 1) 



R 2 



*ia(Ci,C2)=S, a (Ci)^a(C 2 ), 



where £ for principal series of the unitary irreducible representations has the form 

£ — h ip, < p < 00 

2 



(2.78) 
(2.79) 

(2.80) 
(2.81) 



The second order operator Q of eq. (2.21): 

Q = aKl + b(K x K 2 + K 2 K X ) + cK\ + d(K x L 3 + L 3 K ± ) + e(K 2 L 3 + L 3 K 2 ) + jL\, (2.82) 
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can be used to classify all coordinate systems on H 2 . The classification of the operators Q can 
be reduced to a classification of the normal forms of the elements of the Jordan algebra jo(2,l) 
[35] . There are 9 inequivalent forms, in one to one correspondance with the 9 existing separable 
coordinate systems [15, 42, 43]. All the coordinate systems are orthogonal ones. 

The normal forms of the operator Q and the corresponding coordinates are given in Table 
4- Cases I, II, and III correspond to subgroup type coordinates. The corresponding subgroups 
are 0(2), 0(1, 1) and E(l), respectively. The 0(1, 1) subgroup in the equidistant coordinates 
acts in the 01 plane. We could also have chosen the 02 plane (i.e. permuted u\ and u 2 ). 

The elliptic and hyperbolic coordinates of cases IV and V are given in algebraic form. 
Equivalently, they can be expressed e.g. in terms of Jacobi elliptic functions. This makes it 
possible to express the coordinates in the ambient space directly, rather than their squares. 
Indeed, if we put 

Qi — ai — (ai - a 3 )dn 2 (a, k) , g 2 — a± - (ai - a 2 )sn 2 (/3, k') , (2.83) 

and 

e = ^^, k >2 = ^^, k 2 + k ,2 = i, 

ai - a 3 ai - a 3 

into the expressions 

u 2 = r2 (Qi - Q 3)(fe - ag) u 2 = r2 (Qi -fl2)(g2 -a 2 ) u 2 = r2 (Qi - ai)(ai 
(ai - a 3 )(a 2 - a 3 ) ' 1 (a 2 - a 3 )(a 1 - a 2 ) ' 2 (ai - a 2 )(a 1 

we obtain the elliptic coordinates in Jacobi form 

Uq = Rsn(a, fc)dn(/3, k'), u± — i Rcn(a, k)cn((3, k'), u 2 — iRdn(a,k)sn((3,k'), (2.86) 

a G iiK', iK" +2K), (3 e [0, AK'). 



(2.84) 



Q2 



- a 3 ) 



, (2.85) 



3 Contractions of the Lie algebra and Casimir operator 



3.1 Contractions from o(3) to e(2) 

We shall use i?" 1 as the contraction parameter. To realize the contraction explicitly, let us 
introduce homogeneous or Beltrami coordinates on the sphere, putting 

u„ 



X /, 



«3 v /l-( M 2 + w 2 )/jR2 ' 



A* = 1,2. 



(3.1) 



Geometrically (xi,x 2 ) correspond to a projection from the center of the sphere to a tangent 
plane at the North pole. In this parametrization the metric tensor has the following form 



9nv 



r 2 /R 2 



1 



[ii/ 



r 2 /R 2 



(3.2) 



18 



Table 4: Orthogonal systems of coordinate on two-dimensional hyperboloid H 2 



Coordinate systems 


Coordinate 


Limiting 


Limiting 


and integrals of motion 






systems on E 2 


systems on 


I. Pseudo-spherical 


u = 


R cosh r 


Polar 


Cartesian 


t > 0, </? E [0, 2tt) 


Ui = 


R sinh t cos ip 






VS — ^3 


u 2 = 


R sinh r sin cp 






II. Equidistant 


Uq = 


R cosh T\ cosh r 2 


Cartesian 


Polar 




Ui = 


R cosh Ti sinh r 2 






Qeq = Kf 


u 2 = 


R sinh Ti 






III. Horocyclic 


Uq = 




Cartesian 


Rectangular 


y > U, x G IK 

= (Ki + L 3 ) 2 

-t-l \ -L iJ / 


U\ = 

u 2 = 


R(x 2 + y 2 - l)/2y 
Rx/y 




Cartesian coordi- 
nates rotated by 5 
(nonorthogonal) 


IV. Elliptic 


ul = 


J^2 (pi— «3)(P2— 03) 

fai— 03^02— 03^ 


Elliptic 


Elliptic I,II,III 


0-3 < a 2 < P2 < d\ < Pi 


u\ = 


(Pi-a2)(P2-a2) 
fai— a,2)( 0,2-0,3) 


Parabolic 


Cartesian 


Q E = L\ + sinh 2 fK\ 


u\ = 


j£2 (pi-ai)(a 1 -p 2 ) 
(0,1—0,2) (0,1—0,3) 


Cartesian 




V. Hyperbolic 


r.,2 

U Q = 


j^2 (pi—a2)(a2—p2) 
(0,1-0,2) (0,2-0,3) 


Cartesian 


Elliptic II 


P2 < 03 < a 2 < cli < pi 




j^2 (pi-a3)(a3-p2) 
(ai— 03)^02—03) 




Parabolic I 


ts2 -2 t2 

Qh — K 1 — sin aLg 


u\ = 


j^2 (pi-ai)(ai-p2) 
(ai— a2)(ai— 03) 






VI. Semi-hyperbolic 


Uo+ui 
R 2 


— Cl 4- // 2 V1 4- ;/ 2 ^ 

— ^1 "I" A*aA-L "I" 


Parabolic 


Cartesian 


Pl,2 > 


u 2 -u\ 
R 2 


= (1 + pip 2 ) 


Cartesian 




Q SH = -{K 1 ,L 3 } 


u 2 = 


Ry/HlH2 






VII. Elliptic-parabolic 

flGE,i?G (-7r/2,7r/2) 
Qpp = (Kx + L 3 ) 2 + K\ 


Uq = 
Ui = 

u 2 = 


j^cosh a + cos v 
2 cosh a cos j) 

^sinh a — sin d 
2 cosh a cos $ 

i?tan-$tanha 


Parabolic 


Hyperbolic II 


VIII. Hyperbolic-parabolic 


U = 


j^cosh 2 b + cos 2 1? 

2 sinh 6 sin $ 
^sinh b — sin 1? 

2 sinh 6 sin $ 
R cot i? coth 6 


Cartesian 


Hyperbolic III 


6 > o,i? e (o,tt) 

Qpp = (K, + L 3 ) 2 - ki 


Ui = 

u 2 = 






IX. Semicircular-parabolic 


Uq = 
Ui = 




Cartesian 


does not correspond 
to a separable 




u 2 = 


2 r-2 




coordinate system 


HK2, L 3 } 
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The Laplace-Beltrami operator corresponds to 



1 + 



R 2 



d 2 d 1 / d 

dxj + R^dxJ + R 2 \"dx,. 



7TJ + 7T 2 + 



T 2 

2 , ^3 



P 2 



where 



7T,, 



9 (9 



P 2 dxj ' 



- x 2 



_d_ 

dx\ 



Using the connection between operators 7r M and the generators of the 0(3) group 

Li L 2 
~R =n2 R 
we obtain the following commutation relations 



7Ti, L 3 = -(xin 2 - x 2 n 1 ) 



[L 3 , 7Ti] = 7T 2 , [L 3 , 7T 2 ] = -7Ti, [7Ti, 7T 2 ] = -^|, 



(3.3) 



(3.4) 



(3.5) 



so that for P — > 00 the o(3) algebra contracts to the e(2) one. Moreover the momenta ir^ 
contract to P M = d/dx^, (/x = 1,2) and the o(3) Laplace-Beltrami operator (3.3) contracts to 
the e(2) one: 



k L_B 



< + vr 2 z + 



I 2 
R 2 



A = (P 2 + P 2 ). 



(3.6) 



3.2 Contractions from o(2, 1) to e(2) 

As in section 3.1, let us introduce the Beltrami coordinates on the hyperboloid H 2 putting 

= . // 1.2. (3.7) 



x =R^t = R - ^ 

u o ^JR 2 + uj + 



The 0(2, 1) generators can be expressed as: 
Pi d 



= 7T 2 

R 



dxo R 2 



x 2 1 d d 
hi 2 - 



dxi 



' dx 2 
L 3 



~ll= 7Tl=Pl - R 2 



Xi 



d 
dxi 



+ x 2 



_d_ 

dx 2 



Xi7T 2 - X 2 7Ti. 



The commutation relations of the o(2,l) algebra (2.78) in terms of the new operators take the 
form 

[VTI, 7T 2 ] = - — , [L 3 , 7Ti] = 7T 2 , [tT 2 , L 3 ] = 7t\ , (3.8) 

so, that for R — > oo the o(2,l) algebra contracts to e(2) and the momenta 7f M to P M = d/dx^. 
The o(2,l) Laplace-Beltrami operator (2.2) contracts to the e(2) one: 



Alb = 7Ti + 7T2 



r 2 

^3 . A / r>2 



P 2 



A = (P 1 2 + P 2 2 ). 



(3.9) 
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3.3 Contractions from o(2, 1) to e(l, 1) 
Let us introduce Beltrami coordinates on hyperboloid H 2 



u 



yo = R— = R— r 

M 2 ^ul + uj-R 2 
The 0(2, 1) generators can be expressed as 



, Vi = R 



R 



Ui 



u 2 jul + ul-R* 



(3.10) 



~R =7ri 



JL _ y± ( JL _|_ JL 

dy # 2 V °%> 1 dy li 



R~ n<2 dy x R 2 \dyo Vl dyi 



-K 2 = K = y fi 2 + yi7fi- (3.11) 
The commutators of the o(2, 1) algebra (2.78) in the new operators (fri, n 2 , K) take the form 

K 



-VT 2 , 



7Ti. 



(3.12) 



so, that for i? — > oo the o(2, 1) algebra contracts to the e(l, 1) one. The o(2, 1) Laplace-Beltrami 
operator contracts to the e(l, 1) one: 



,2 ,2 K 2 d 2 d 2 
Alb = TTi - tt? H -* 



(3.13) 



2 ' i? 2 ' dyi dyV 

and eq. (2.80) transforms for large t ~ pR to the one-dimensional Klein-Gordan equation. 

(3.14) 



d 2 ib d 2 ib 9 



dyg <% 2 

4 Contraction for systems of coordinates 
4.1 Contractions and coordinate systems on Si 

1. Spherical coordinates on S 2 to polar on E 2 . 

We consider the spherical coordinate (2.64) with the parametr a\ = a 2 and put 



tan6> 



R 



In the contraction limit R — > 00, — > we have 



xi = it > x = r cos ip, x 2 = R > y — r sm ip 

u 3 u 3 
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and 



Qs = L 



2. Spherical coordinate on S 2 to Cartesian on E 2 . 

We choose the spherical coordinate (2.64) with a 2 = a 3 . Putting 

cos 6' ~ — ~ 0, cos ip' ~ ^- ~ 0, 

fx it 

and taking the limit R — > oo and — > |, <// — > |, we obtain 

1 L 2 

^^■5 = ^ = ^ P i 2 ~ @c 

and 

cot 0' „ . 

x\ — R > x x 2 = R cot if — > y. 

sin <// 

It is easy to see that for the case a± = a 3 the corresponding spherical system of coordinates 
(2.64) contracts to Cartesian coordinates on E 2 for R — > oo. 

3. Elliptic coordinates on S 2 to elliptic coordinates on E 2 . 

We take Q in its general form, equivalent to 

Q E = Ll-( a -^)L\. (4.1) 



We put 



a 3 - a 2 



R 2 D 2 



a 3 — ai a 2 — ai 
and in the limit R 2 ~ a 3 — > oo obtain 



(4.2) 



Q £ = L 2 - ^ L 2 - L 2 - D 2 P| ~ J*. (4.3) 
\a 3 — a 2 / R z 

For the coordinates we put 

Pi — ai + (a 2 - ai) cos 2 77, p 2 = «i + (a 2 - ai) cosh 2 £, (4.4) 

and for i? 2 ~ a 3 — > 00, using eq. (4.2), we obtain eq.(2.41), i.e. elliptic coordinates on the plane 
E 2 . 

4. Elliptic coordinates on S 2 to Cartesian coordinates on E 2 . 
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We start from the coordinates (2.72) but change the ordering of the parameters 
corresponds the interchange of coordinates u 3 <-> u 2 , i.e. put 

oi < Pi < a 3 < p 2 < a 2 , 

and choose — a\ = a 2 — = a. Than we have 

Q E = a(Ll-Ll). 

Introducing the new coordinates by 

a 3 ~ Pi t P2 - % 

- 6 - 6, 



we can rewrite the (2.72) in the form 
R 2 



R 2 



< = -y(i - 6)(i + 6), < = — (i + 6)(1 - 6), «2 = ^66- 

Using eq. (3.1) we have for Beltrami coordinates 



x 2 = ^(i - 6)(i + 6) _a _ D2 (i + 6)(i - 6) 



266 



= R A 



266 



^From equation (4.9) we obtain 
6,i 



R 2 + x\ + x\ 
Taking now the limit R — > oo we have 

6^1- 



i? 2 " l/, >; 



1/2 2 2 
1 2 



2R 2 



x 



1 - 



R 2 ~ ' ~ R 2 ' 
and hence x\ and rr 2 of eq. (4.9) go into Cartesian coordinates: 

xi — > rr rc 2 — > y. 

For the integral of motion in the limit i? 2 ~ a — > oo we have 

1 



aR 2 



Q E = {K-K)^Pt-P2=Qc- 



5. Elliptic coordinates on £2 to parabolic coordinates on E 2 . 
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We take the operator (2.72) with a\ < p\ < a 2 < p 2 < ^3 and choose the parameter 
03 — a 2 = a 2 — ax = a. We must first "undo" the diagonalization (2.26) by a rotation through 
7r/4. The operator (2.63) transforms into 



sQe = — ^(LiLz + L3L1) = (L 3 7r 2 + 7r 2 L 3 ), 
art ti 



(4.14) 



with the correct limit (2.39) for R — > 00. The coordinates (2.77) on ^ are rotated into 
Mi 



-^(snadn/? + dnasn/5), «2 = -Rcnacn/?, w 3 = -^(dnasn/? — snadn/5), (4.15) 
V 2 y2 



with modulus k = k' = 1/ \[2 for all Jacobi elliptic functions. 
/,From eq (4.15) we obtain 



sna 



v^dn/? = 



1 

71 
1 

72 



1/2' 



^y /2 f 1 + ^y /2 



R) 



Eq. (4.16) suggest the limiting procedure. Indeed we put 



2 2 
sna = -l + ^, v^dn/3 = l + |^. 



In the limit R — ■> 00 we obtain 

X! — > X = 

i.e. the parabolic coordinates (2.39). 



■u 2 - 



x 2 -> y = «u, 



(4.16) 



(4.17) 



(4.18) 



4.2 Contractions of coordinate systems from #2 to £3 



1. Pseudo-spherical Coordinates on H 2 to Polar Coordinates on E 2 

In the limit R — > 00, r — >■ putting tanhr ~ r/i? we have: 



r2 

-^3 



r 2 

-^3 



and for Beltrami coordinates (3.7) we obtain: 



xi = -R- 



U 



x = r cos v?, x 2 = R 



U2 
U 



y — r sm 93. 



2. Equidistant Coordinates on if 2 to Cartesian on E 2 
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For Beltrami coordinates (3.7) we have: 

xi = Ptanhr 2 , x 2 = Ptanhri/ coshr 2 . (4.19) 

Taking the limit P — > oo, Ti,t 2 — > and putting sinhrx ~ y/P, sinhr 2 ~ x/R in (4.19) we 
obtain x± — > x, x 2 — > y and 



P 2 



= < -> p/ ~ g c . 



3. Horocyclic Coordinates on P 2 to Cartesian on P 2 

For variables x, y we obtain: 

u 2 _ P 
x = , y 



M — Mx M — 

In the limit P — > oo we get: x — > y/P, y — > 1 + x/P and Beltrami coordinates go into 
Cartesian ones 

n i 2 + § 2 -l 2xP 
an = P -> X, x 2 = -> y. 

x 2 + y z + 1 x 2 + y z + 1 

For integral of motion we have: 

^ 2 - = - 2 + ^--{ 7 r 2 ,L 3 }-P 2 ~g c . 

4. Elliptic Coordinates on P 2 to Elliptic Coordinates on P 2 

We put 

P 2 D 2 



(4.20) 



a 2 — a 3 oi — a 2 
and in the limit P 2 ~ (—03) — > 00 obtain: 

^ = L l + - ^ 2 + ^ 2 Pi ~ Qe, 

where 2D is the focal distance. Writting the coordinates as 

Pi = a,i + (ax - a 2 ) sinh 2 £, p 2 = a 2 + (ai - a 2 ) cos 2 77 

and using eq.(4.20) in the limit P 2 ~ (—03) — > 00 we get the ordinary elliptic coordinates on 
P 2 plane [15, 10] 

5. Elliptic coordinates on P 2 to Cartesian on E 2 . 
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We make a special choice of the parameters af a\ — a 2 = a 2 — a 3 and determine new variables 
6,2 by the formula 



g h2 -a 2 _ul + uj 

6,2 - - — rrr^- ± 



a\ — a 2 



2R 2 



' 2 i 2 \ ^ 
«0 + U 2 \ 



\ \ 2R 2 



R 2 ' 



(4.21) 



Considering the limit R — > oo we obtain: 6 ~ 1 + 22/ 2 /i? 2 , 6 

~ x 2 /R 2 and the Beltrami 

coordinate (3.7) take the Cartesian form 



Xi 



R^ 

u 



R\ 



266 



(6 + 1) (6 + 1) 



X, x 2 



R^ 

U 



(6-i)(i-6) 
^(6 + i)(6 + i) 



The operator Q E goes to Cartesian one 



Qe _ L\ 2 
R 2 ~ R 2 + 711 



6. Elliptic coordinates on if 2 to parabolic on i? 2 . 
We start from the rotated elliptic coordinates 




/ cosh / sinh / 
sinh / cosh / 
V 1 




/ Mo cosh f + ui sinh / 
Uq sinh / + ui cosh / 
V u 2 



(4.22) 



where sinh 2 / = (ai — a 2 )/(a 2 — a 3 ). We choose a 2 — a 3 = ai — a 2 = a. Then for rotated elliptic 
coordiantes (4.22) we get 

Mo — —i=(snadnf3 + i\/2cnacnf3) , wi = -^(icnacn/? + v^sncudn/?) , w 2 = iRdnasnfi, (4.23) 
V 2 V 2 

with modulus k — k' — 1/ \[2 for all Jacobi elliptic function. The integral of motion transforms 
into 

Q*v = 3L 2 -V^i^Ls + L 3 #i), (4.24) 
with the correct limit to (2.36). From eq. (4.23) we obtain 



cna 



1 + 



Ur 



+ zrh + 



2\ V W2 i?/ 2i? 2 ^2\V Ry/2 R 



1 - 



Ur 



2i? 2 ' 



cn/3 



1 + 



2\\^RV2 Rj^2R 2 ^2\\ RV2^ R 



Un 



2R 2 ' 
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and therefore for large R we have 

1 u 2 Q n 1 v 2 

—tcna ~ 1 = — , cnp ~ 1 H = — . 

2^2 i? 

In the limit R — > oo we obtain 

2 2 

■u - V 

xi -> x = — - — , x 2 -> y = «i>, 
i.e. the parabolic coordinates (2.39). 

4.3 Contractions of coordinate systems from H2 to £^1 

1. Equidistant coordinates on H2 to pseudo-spherical ones on E\ t \ plane. 
For Beltrami coordinates (3.10) we have: 

y = -Rcothri coshr 2 , y\ = .Rcothri sinhr 2 . (4-25) 

Taking the limit R ^ 00, Ti^^f + Jj and putting 

coth n = tanh ^ ~ (4.26) 

we obtain 

y — > t = r cosh r 2 , yi — >■ x = r sinh r 2 , (4.27) 
where < r < 00, — 00 < r 2 < 00. For the integral of motion we get 

Qeq = K\ ^Q s = L\. (4.28) 

2. Pseudo-spherical coordinates on H 2 to Cartesian coordinates on E± t i 
For coordinates (3.10) we have 

cothr „ . , , 

y = R , yi = Rcot<p. 4.29 

COS (f 

Taking the limit R^oo,T^in/2,(p^^ and putting 

t x 

cothr~— , coty?~— , (4.30) 
R R 

we see that Beltrami coordinates go into Cartesian ones 

2/o -> t, yi -»• x. (4.31) 

For the integral of motion we obtain 

t = «c (4.32) 



27 



5 Contraction of basis functions on S2 and H2 

5.1 Contraction of spherical basis and interbasis expansions 

5.1.1. Spherical basis on 5*2 to polar on E 2 . 

We start from the standard spherical functions Yi m (6,<f)) as basis functions of irreducible 
representations of the group 0(3) (see e.g. ref. [37]) 



m+ |m I 

Y lm (9,<f>) = (-1) — 



2/ + l(/+|m|)! 
2 (l-\m\)\ 



1/2 (sin 



2H|m|! 



2 Fi [-1 + |m|, I + \m\ + 1; \m\ + 1; sin 2 - —= (5.1) 



2) v/27 



In the contraction limit R — > oo we put 



tan#~#~-j-, I ~ kR. (5.2) 
R 



Using the asymptotic formulas 



and formula 



/ r 2 \ / k 2 r 2 \ 

l-kR,kR;\m\ + l; — \ = F 1 (\m\ + 1; — J (5.3) 

\im^±4 = (5.4) 



we obtain 



1 mi „imd> 

lim -= Y lm {9,<j>) = (-l)^-VkJ H (kr)-= (5.6) 

ft^oo VR V27T 

The result (5.6) is not new [37]. The point is that this asymptotic formula is obtained very 
naturally in the context of group contractions applied to the separation of variables. 

5.1.2. Spherical basis on S 2 to Cartesian on E 2 

We start from the coordinates (#', 0') in eq. (2.64), but drop the primes, and write the 
corresponding spherical functions as 

Y lm {e,<t>) = ^n e ^ (sin 0)M 
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x < 



l + m 
(-1) 2 



r ( i±m±l )r( i=m±l ) 
r(( i + m + 2 )r( i-m + 2 ) 



2^1 



/-to Z+m+1 . 1. rri q2 , 
2 ' 2 ' 2 ' 



(5.7) 



i+m-i r r(' + ™ +2 )r('~^ +2 ) 

(( i±m±l )r( tm±l ) 



. tT'"-i r 1 

"I) 2 [f 



* 2 cos 02^1 (-^f^,^f±^;|;cos 2 ^ 



for / + m even and odd, respectively, we now put 

I ~ ZcR, m ~ k 2 R, 

and 

sin# — > 1, cos6> — > ^, cos0 



7T 

2' 



2/_ 



(5.8) 



(5.9) 



The 2 -Pi hypergeometric functions simplify to o-Pi ones, the T functions also simplify and the 
final result is that under the contraction we have 

rJ2_2 - 



lim (-1) 

H— >oo 



f-j-m 



7T 



e 



k -i(fc 1 x) F 1 (|; ^ 

COS fciX 



• < 



k —i sin fciX 



(5.10) 



with k\ + A; 2 2 = k 2 and for Z + m even and odd, respectively. The parity properties of Y\ 



Irn 



under the exchange 9 — > 7r — have lead to the appearance of cos/cix and sin/cix in eq.(5.10), 
instead of the usual Cartesian coordinate solution expi(kix + k 2 y). 

Finally note that, the factor in formula (5.10) are connected with the contraction of 
Kronecker symbols to delta function 

^^^(k-k') = ^5(k l -k' l ) 



5.1.3. Contraction in interbasis expansions. 

Let us now to consider the contraction R — > oo for the interbasis expansion (2.69). The 
contraction of basis functions were presented in the formulas (5.6) and (5.10). In order to obtain 
the corresponding limit we need the asymptotic behavior of the "little" Wigner <i-function for 
a large R. It is easy to see that the expression of rf-function for the angle n/2 in terms of 
hypergeometric functions 2 Fi (see for example [37]) is not applicable for the contraction limit 
when £ — > oo and m — > oo simultaneously. To make this contraction we use an integral 
representation for the function d l mi (P) [37] 

i 

(5.11) 



Cm!*/*) = 



■m2—m\ 



2tt 



(l + m 2 )\(l -m 2 )\ 
(l + mi)\(l - mi)! 
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X 



2tt 



e* 2 cos 



— h le 2 sm- 
2 2 



| le- 1 * cos | + ie* a sin | e im2V efy>, 



which for the particular case of /? = 7r/2 can be presented in the following form 

' (I + m 2 )\{l - m 2 )\ 1 



di 



, 7T . 



I-m t 2' 

7T 



(Z + 77li)!(Z-77li)! 



sin ay~ mi (cos a)' +mi e 2im2Q cfa (5.12) 



Using now the formulas [44] 

cos(2na) = T„(cos2o;), sin(2no;) = sin 2a ■ C/ n _i(cos 2a), 

where T;(x) and C/j(x) are Tchebyshev polynomials of the first and second kind. After integrat- 
ing over a, we obtain a representation of the Wigner <i-function for angles tt/2 in terms of the 
hypergeometrical function 3-F 2 (l) 



dU,mi(|) = ^^V / ( / + ^)!(/-m 2 )! 



(5.13) 



x < 



r(i±=i+i)r(^Ti) 



3-^2 



1+mi+l 
"777,2, r«2, g 



i / + 1 



(/ + mi) — even, 



2ii 



rfi±jii+i)r(^i+^ "i - 



[ (i+i) \ r(i±^±i)r(i=^H) 
For large i? we put 



3-^2 



-7712 + 1, 7712 + 1,^ + 1 

|,/ + 2 



1 , (I + m 1 ) - odd. 



l~kR, mt-hR, & X ~1L, 2 ~-|, 



(5.14) 



where k 2 = k\ + fcf. Using the asimptotic formulas for 3-F 2 (l) function (5.3) and T function 
(5.4), we get 



/—mi , , 7T / / 



x < 



(£); 2 F 1 (-m 2 ,m 2 ;|;^), = ( _ 

[ -im 2 (ff) 1 2 Fi (-m 2 + l,m 2 + 1; §; ^) , 
with cosy? — ki/k and for (/ + mi) even or odd respectively. 



1) : 



Tlko 



cosm 2 y2, 
k isinm 2 (/9, 



(5.15) 
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l~ l m l I 

Multiplying now the interbasis expansion (2.69) by the factor (— 1) 2 ; and taking the 
contraction limit R — > 00 we obtain (6 = 6 2 , m = m 2 ) 

^{^I'l^l-Z}^" - ' (5 ' 16) 

or in exponential form 

00 

e ikrcos(6-<p) = (i) m J m (kr) e im{e -^ . (5.17) 

m=— 00 

The inverse expansion is 

J m {kr)e ime = ± — '— e im V -ikr cos (e- v ) d(p ^ ^ 

2ir Jo 

For 6 = the two last formulas are equivalent to the well know formulas in the theory of Bessel 
functions [44], namely expansions of plane waves in terms of cylindrical ones and vice versa. 

5.2 Solutions of the Lame equation 

Let us consider eq. (2.60) on the sphere S 2 and separate variables in the elliptic coordinates 
(2.72). We obtain two ordinary differential equation of the form 

d 2 tp | 1 f 1 1 1 1 # 1 f A - 1(1 + l)p 

dp 2 2~\p-ai p-a 2 p - a 3 j dp A{(p-a 1 )(p-a 2 )(p-a 3 ) 

or equivalently 
where 

p {p) = (P- a i)(p ~ a 2){p ~ a 3 ). 

Eq.(5.19) is the Lame equation in algebraic form. It is a Fuchsian type equation with 4 regular 
singularities (at ai, a 2 , a 3 and 00) [38, 39, 40, 45, 46, 47]. 

Its general solution can be represented by a series expansion about any one of the singular 
points at as 

00 

1>{p) = (p-a i r/ 2 (p-a 2 r^ 2 (p-a 3 r^Y, b{ t k \p-a k ) t , (5.21) 

t=o 

where we have 

«i(«i- 1 ) = > J = 1,2,3 



i> = (5.19) 
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and can choose k equal to 1,2, or 3. 

Substituting (5.21) into the Lame equation (5.19) we obtain a three term recursion relation 
for b\ 

Pt k) b?+i + hi k) + A - Z(Z + l)a k \bf ] + (2t + a-l-2)(2t + a + l- l)bf\ = (5.22) 

with 

a = ai + a 2 + a 3 , a ik — - a k , 6_i = 

A W = 4(a i -a fc )(a J --a fc )(t+l)(t + a fc + l/2) (i, j, fc cyclic) 

7t (fc) = - a k )(2t + a k + aj) 2 - (aj - a k )(2t + a k + ai) 2 . (5.23) 

The expansion (5.21) represents a Lame function. Since we are interested in representations of 
0(3), the sum in ip(p) must be a polynomial of order N, i.e. we must have 

b N ^0, b N+1 = b N+2 = ■■■ = () (5.24) 

for some N. The condition for this is that we have 

/ = 2N + a (5.25) 

and we obtain a secular equation for the eigenvalues A, i.e. the separation constant in elliptic 
coordinates, by requiring that the determinant of the homogeneous linear system (3.14) for 
{b , bi, • ■ ■ , &jv} should vanish. Since iV and I must be integers, eq. (5.25) implies that a and / 
must have the same parity. 

Numerous further properties of the Lame polynomials, in the context of representations of 
the group 0(3), in the 0(3)dD2 basis, were established e.g. in Ref. [38, 39, 40]. 

Here let us just represent the basis functions as 

^(Pi,P2)=AM(Pi)C(P2), (5.26) 

where Af£ is some normalization constant. The labels p, q take values ±1 and identify represen- 
tations of D 2 . For each value of / the values of p, q and A label (21 + 1) different states. Since a 
given representations (p, q) of D 2 can figure more than once in the reduction of a representation 
of 0(3) corresponding to a given /, we are faced with a "missing label problem", resolved by 
the quantum number A, i.e. the operator Q of eq. (2.63). 

The expansions that we shall use for the Lame polynomials in (5.26) are as in eq. (5.21), 
but the summation over t is from t — to t — N. 
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5.3 Elliptic basis on S2 to Cartesian basis on E2 

We choose elliptic coordinates on S 2 as in eq. (2.72), but with a± < a 3 < a 2 , as in eq. (4.5). We 
write the basis functions as in eq. (5.26) with 

N 

Mpi) = (pi - «i) qi/2 (pi - a 2 r /2 ( Pl - a 3 r /2 5> (1) (pi _ ai y 

t=0 
N 

MP2) = (P2 - a i r/\p 2 - a 2 r/\p 2 - a,r^Y, h f\p2 - 02)' (5.27) 

t=o 

as in eq. (5.21). The coefficients (j = 1, 2) satisfy the recursion relation (5.22) and we have 
N = (I — a)/2. We use the coordinates £1 and £ 2 introduced in eq. (4.7) (for a = a 3 — a\ = 
a 2 — a 3)- Eq- (5.28) reduces to 

, 3 N 



Mb) = (-l)^af(l-6)^(l + ei)^^E^( 1 -6)' 

t=o 

N 

m&) = (-ij^a-^^a+^^E^ 1 -^ ( 5 - 28 ) 



«1 . S3. 

t=0 



with C t = afbt, C\ = (— a)*6 t . The recursion relations (5.22) now imply 

8(t + l)(t + e*i + ^Cffi + - 2(2* + ai + a 3 ) 2 - (2t + «i + a 2 ) 2 }C t (1) 

+(2t + a-l-2)(2t + a + l- 1)C^\ = 
-8(t + l)(t + a 2 + \)Ci% + {/i (2) + 2(2t + a 2 + a 3 ) 2 + (2t + a x + a 2 ) 2 }Cf ) (5.29) 

-(2t + a - Z - 2)(2t + a + I - l)C[% = 

where 

fi {j) = -[\ - djl(l + 1)], j = l,2 (5.30) 

(X 

The contraction limit is taken using eq. (4.7) to relate £i j2 to the Cartesian coordinates on 
E 2 . Taking Z ~ kR we find 



fi^->2R 2 kl ^^-2R 2 k 2 2 , k = y/kl + k 2 2 . (5.31) 



For R — > 00 the recursion relations (5.29) simplify to two term ones that can be solved to 
obtain 
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«7 + 77 = \OCj + - «j + - 



with 

Substituting (5.3) into (5.28) we obtain 

Mti) 



t>\ [a-j + - 



1. 



ao+ao Q2 I { 

(-l)^^—x a \F 1 «! + -: 



lM6) = (-l)^^ Q2 o^i(«2 + ^: 
Using now the formula (5.10) we find the contraction limit: 



^(-l)^af x 



cos /cix cos A; 2 y, 
— ^ cos k\X sin /c 2 |/, 
sin /cirrcos /c 2 ?/, 



«i = 0, Q.2 — 
«1 = 0, tt 2 = 1 

d\ = 1, a 2 — 



k\k2 



^2 sin /cix sin fc 2 ?/, «i = 1, a 2 = 1. 



(5.33) 



(5.34) 



5.4 Elliptic basis on 5*2 to elliptic basis on 

Let us start from the elliptic coordinates (2.72) with a\ < p\ < a 2 < p 2 < 03. We take 
the limit R — > 00, a 3 — > 00 with y/a^/R, a± and a 2 finite. We introduce a constant .D as in 
eq. (4.2). Elliptic coordinates on the plane E 2 are introduced via eq. (4.4), so that the Cartesian 
coordinates (x, y) are expressed in terms of the elliptic ones (£,??) as in eq. (2.41). Let us first 
take the limit in the separated equations (5.19). Going over to the variables (£,??) from (pi,p 2 ) 
we obtain for R — > 00: 



drf 



+ \V 



k 2 D 2 f a 2 + a x 
a 2 - ai 



k 2 D 2 



cos 2i] > ipi = 



(5.35) 



«P*,f. W^ + M _^ cosh2 a^ 2=0 (5.36) 




with 

/i = — , / ~ fci?. 

a 3 

In (5.35) we recognize the standard form of the Mathieu equation, whereas eq. (5.36) is a 
modified Mathieu equation [41]. Thus, in the contraction limit, Lame functions will go over 
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into Mathieu ones. Moreover, periodic solutions of the Lame equation go over into periodic 
solutions of eq. (5.35). 

The contraction limit can also be taken directly in the Lame polynomials, using the expan- 
sion (5.21) (cut off at t — N). The result that we obtain is 

tL^T = (« 2 -«i) f ^f^(cosr ? )-(sin^)-£a(cosr ? ) 2 ', (5.37) 

^^^ = (a 2 -aO^^^(coshO a HsinhO a2 E^(^h0 2 ', (5.38) 

where the expansion coefficients C t satisfy recursion relations obtained from eq. (5.22), namely 
4(t + \){t + 1/2 + ai )C t + {/i - (2t + ai + a 2 ) 2 }C t - k 2 D 2 C t = 0. (5.39) 

5.5 Elliptic basis on S2 to parabolic basis on E2 

Let us consider the contraction limit for the Lame equations (5.19). To do this we use equations 
(2.76) with a 3 — a 2 = a 2 — a\ = a i.e. k — k! — l/y/2, together with eq. (4.17), to obtain 

Pi ~ ai + a ^-1 + , p 2 ~ ai + a ^1 + ^ (5.40) 

The equation (5.19) for p — p\ and p = p 2 in the limit i? — > 00, with / 2 ~ /c 2 _R 2 and A — a 2 l(l + 
1) = pRa, yields the two equations 

+ {k \2 + = Qj ^ + ^2^2 _ ^ = (5 41) 

respectively. 

Thus the Lame equations in the contraction limit go over into the equations (5.41) for 
parabolic cylinder functions [44]. The same is of course true for solutions. The expansion 
(5.21) is not suitable for the contraction limit. In view of eq. (4.17) we need expansions in 
terms of the variables (1 + sna) and (1 — \/2dn(3). This is not hard to do, following for instance 
methods used in Ref.[48] to relate the wave functions of a two-dimensional hydrogen atom, 
calculated in different coordinate systems. The formulas are cumbersome, so we shall not 
present them here. 

5.6 Contractions of basis functions from H2 to E2 and 
5.6.1. Pseudo-spherical basis on H 2 to polar basis on E 2 . 
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The pseudo-spherical eigenf unctions \J/ pm (r, <p) normalized to the Dirac delta- function, have 
the form: 



*p"»(t,vO = I r(i + tp+ H) | -pHjcosh^exp^), (5.42) 



2 



where m = 0, ±1, ±2, .... In the contraction limit R — > oo we put: tanhr ~ r ~ r/i?, p ~ A;i?. 
Rewriting the Legendre function in terms of hypergeometric function as [49] 

P 'i(coshr) = ( . —-7 , , 2-Fi ( o + H+ *P>o + m m ;-smh - 

l P 2 r(i + ip — \m\) |m|!2l m l \2 2 2 

Then using the asymptotic formula for hypergeometrical function 2 F\ and Y function 

lim | T(x + iy) | exp(- | y |) | y \?~ x = y/2n (5.43) 
\y\— >°° 2 

we obtain in the contraction limit R — > oo: 



lim ^ pm (r,<p) = v 7 ^- J| m |(Ar) • -7=. 



5.6.2. Pseudo-spherical basis on H 2 to Cartesian basis on E lj± . 

Taking the Legendre function in eq. (5.42) in terms of two hypergeometric functions [49] 

P m (cQshT) ^(sinhr)- L ^(f ~ lg)£(! ~ gg) 

,_, /3 m + ip 3 m — 3 , 9 \ ^ /l m + ip 1 m — ipl , 9 

_p -, -; -: cosh 2 r + 2-F1 -, -; -; cosh 2 r 

V4 2 '4 2 ' 2' / V4 2 '4 2 ' 2' 

Putting for large R 

p ~ kR, m ~ kiR, cothr ~ — , cot<p ~ — , ; k 2 + k\ = k 2 . 

R R 

Using two asymptotic formulas 

,_, (\ m + ip 1 m — ipl l9 \ „ / 1 ^ 2 £ 2 \ /, x 

& 2Fl (1 — r 2 - 4 — i ; cosh T ) = » Fl U ; - -rj = cos( *'> t) - 

, ^/3 m + ip 3 m — ip 3 l2 \ / 3 & 2 £ 2 \ sin(fcoi) 

hm 2-P1 7 n — , 7 n — ; 7; cosh t)= F 1 \-; ' 



«-«, x V4 2 ' 4 2 ' 2' J u ^ V 2 4 / fc ot 

and formula (5.4) we finally obtain 

lim v^|r(ip)|* pm (r, <p) = & fe °'-**. (5.44) 
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5.7 Contractions for equidistant basis on H2 

5.7.1. Equidistant basis on H 2 to Cartesian basis on E 2 . 

In the equidistant system the normalized eigenfunctions ^ p \(ti, t 2 ) have the form: 



To perform the contraction we write the Legendre function in terms of hypergeometric function 
[49] 

^(-tanhn) = ^ - ^ - J _ | 2 F, (- + a, - + b; -; tanh 2 r,) 

+2tanhn r§-l)r(l-!) 2i ^ G (L T k I tanh2r 0}' 

where a = i(p — A)/2, b = i(p + A)/2. For large i? we put p ~ fci?, A ~ r 2 ~ ri ~ 

y/i? where x,y are the Cartesian coordinates. Then using the asymptotic formulas: 



/111 \ /l y 2 k 2 \ 

km Q 2F 1 f - + a, - + b; -;tanh 2 rij = F 1 I -; ^ J = cos fc 2 y, 

il^ 2Fi (i- a 'i- 6; l tanh2ri ) = ° Fi (|-^r) = ^ sinA;22/ ' 

where kj + k 2 = k 2 , we finally get 



lim ^a(ti,t 2 ) = J-^-exp(ik 1 x + ik 2 y). 



5.7.2. Contraction from equidistant basis on H 2 to polar on Ei^. 
Writting the Legendre function in terms of hypergeometric functions [49] 

i^_ 1/2 (tanhn) = ^(sinhr!)^- {2- iA (cothr 1 ) iA+1 / 2 — T Fi ' Vi 



F\ g - ? - 1 + <A; coth 2 n ) + 2*(coth 

r(2A) F.fi-^^-^l-^coth 2 ^!. 



r(|-i(p-A)' l \A 2 '4 
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Putting for large R: p ~ kR and coshri ~ ^, and using the asymptotic formulas for hyperge- 
ometric functions 



& (i - ! - 1 - * coth2 *) = ^ - a > (I) + " j - (fcr >- 

we obtain 

where (z) is the first kind of Hankel function. 
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